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Abstract. We present here generalization of the recursion method of Haydock et 
al [1] for the calculation of Green matrices (in angular momentum space). Earlier 
approaches concentrated on the diagonal elements, since the focus was on spectral 
densities. However, calculations of configuration averaged response functions or 
neutron scattering cross-sections require the entire Green matrices and self-energy 
matrices obtained from it. This necessitated the generalization of the recursion method 
presented here with examples. 



PACS numbers: 71.20.Be, 71.23.-k, 72.15.-v 
1. Introduction 

The augmented space recursion carried out in a minimal basis set representation of the 
tight-binding linear muffin-tin orbitals method (TB-LMTO-ASR) has been proposed 
earlier by us [2, 3] as technique for the incorporation of the effects of configuration 
fluctuations for random substitutionally disordered alloys. This can be achieved 
without the usual problems of violation of the Herglotz analytic properties [4] of the 
approximated configuration averaged Green functions for the Schrodinger equation for 
these random alloys. Although our initial focus was with configuration averages of 
density of states and spectral functions, recently we have proposed using the TB-LMTO- 
ASR for the study of configuration averaged optical conductivities [5] or coherent and 
incoherent neutron scattering cross-sections [6] . These calculations require the full Green 
matrices in angular momentum space and not only their diagonal elements. We propose 
here the use of a generalization of the recursion method of Haydock et al [1]. The block 
recursion technique had been introduced earlier by Godin and Haydock [7, 8] in the 
very different context for obtaining the scattering S-matrix for finite scatterers attached 
to perfect leads. We shall borrow their ideas and set up a block recursion in angular 
momentum space (rather than the lead space, as in Godin and Haydock's work) in order 
to obtain the Green matrices (in angular momentum space) directly. 
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2. Methodology 

2.1. The TB-LMTO Hamiltonian in augmented space and the recursion method 

The augmented space recursion based on the tight-binding linear muffin-tin orbitals 
method (TB-LMTO- ASR) has been described thoroughly in a series of articles [9]- [17]. 
We shall introduce the salient features of the ASR which will be required by us in our 
subsequent discussions. 

We shall start from a first principle tight-binding linear muffin-tin orbitals (TB-LMTO) 
[18, 19] in the most-localized representation (a representation). This is necessary, 
because the subsequent recursion requires a sparse representation of the Hamiltonian. 
In this representation, the second order alloy Hamiltonian is given by, 

H (2) = E u + h - hoh 

where 

h = E (Cfl - E vfl ) P^ + EE A R 2 S RR' A r" T RRi 
R R R' 

o = J2o r v r (i) 

R 

C R , E u r , Ar and o R are diagonal matrices in angular momentum space : 

Cr = C RL 5 LL , E uR = E vRL 5 LL , A R = A RL 5 LL , and o R = o RL 5 LL , 

and Srr> = Sr L ,r>l> is a matrix of rank L max . V R = \R){R\ and T RR > = \R)(R'\ 
are projection and transfer operators in the Hilbert space Ti spanned by the tight- 
binding basis {l-R)}. Here, R refers to the position of atoms in the solid and L is a 
composite label {£,m,m s } for the angular momentum quantum numbers. C, A and 
o's are potential parameters of the TB-LMTO method; these are diagonal matrices in 
the angular momentum indices, o _1 has dimension of energy and E^'s are the energy 
windows about which the muffin-tin orbitals are linearized. 

For a disordered binary alloy we may write : 



C RL =C£n R + Cf (1 - n R ) 

^1= (Ai) 1/2 n R+ (A?) 1/2 (l-n R ) 

o RL =o J £n R + of (l-n R ) (2) 

Here {n R } are the random site-occupation variables which take values 1 and depending 
upon whether the muffin-tin labeled by R is occupied by A or 5-type of atom. The 
atom sitting at {R} can either be of type A (n R = 1) with probability x or B (n R = 0) 
with probability y. The augmented space formalism (ASF) now introduces the space of 
configurations of the set of binary random variables {n R } : $. 
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In the absence of short-ranged order, each random variable n R has associated with it 
an operator whose spectral density is its probability density : 

p(n R ) = x5(n R - 1) + y6(n R ) 

= - - lim Im(T* | [(n R + i5)I - M^ 1 | U) (3) 

where Mr is an operator whose eigenvalues 1, correspond to the observed values of 
n R and whose corresponding eigenvectors { | \0 r )} span a configuration space (f> R of 
rank 2. We may change the basis to {| | l R )} : 

I til) = y/x \0 R ) + ^y\l R ) 

I W = Vv \°r) - \i R ) 

and in this new basis the operator is 

n R - M R = xV R + yV l R + ^ (Tr + Tr) (4) 

These two vectors span the space (p R . The full configuration space $ = fll 4>R is 
then spanned by vectors of the form | TTITI • • •)• These configurations may be labeled 
by the sequence of sites {C} at which we have a j. For example, for the state just 
quoted {C} = |{3, 5, . . .}). This sequence is called the cardinality sequence. If we define 
the configuration | ft • • • T • • •) as t ne average or reference configuration, then the 
cardinality sequence of the reference configuration is the null sequence {0}. 

The augmented space theorem [9] states that 

«A({n i? })» = <{0}|A|{0}> (5) 

where 

A({M R }) = J ...JA{{X R }) n^P(A«) 
P(A#) is the spectral density of the self-adjoint operator M^. 

Applying this idea we may obtain exact expressions for configuration averages of 
Green matrices (in angular momentum space) both in the real and reciprocal space 
representations 

< G(R, R, z) > = (i? O {0}|(-2l - H (2) ) _1 |i? O {0}> (6) 

«G(k,z)»= (k®{0}|(^-H( 2 )) _1 |k®{0}) (7) 

where G and H^ 2 ) are operators which are matrices in angular momentum space, and 
the augmented k-space basis |k, L <g> {0}) has the form 

(1/y/N) ]T exp(-ik • R)\R, L ® {0}) 

R 



The augmented space Hamiltonian H^ 2 ^ is constructed from the TB-LMTO Hamiltonian 
jj( 2 ) by replacing each random variable n R by operators M^. It is an operator in the 
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augmented space \& = 7i <8> $. The ASF maps a disordered Hamiltonian described in a 
Hilbert space 7i onto an ordered Hamiltonian in an enlarged space where the space 
^ is constructed as the outer product of the space 7i and configuration space $ of the 
random variables of the disordered Hamiltonian. The configuration space $ is of rank 
2^ if there are N muffin-tin spheres in the system. Another way of looking at is to 
note that it is the collection of all possible Hamiltonians for all possible configurations 
of the system. 

A little mathematics yields the following : 

<G(z)> = (1| (A + B + F-S + (J + S) 6 (J + S)) _1 |l). (8) 
for real space calculations : 

W = WK^w lmm + J^W lR 8 {R}} - 

or for reciprocal space calculations : 

I 1 ) = [A (A-i)]^ |k W) + [A(A-i)]i>« |k {R}) - 



where 



K = E{K((^I-C)A- 1 )/A(A- 1 )} O k 



R 

and K may be A, B or F, while the operators : 

O a = V r ®1 O b = V r ®V r F = V R ®{T£ i + T£} 
Moreover, J = J a + Jb + Jf and 6 = 6^ + 6 B + 6 F where : 



J K = E{K((C-E,)A- 1 )/A(A- 1 )} O k 

R 

6 K = E{K(oA) A(A- 1 )} O k (9) 

R 

For any diagonal (in angular momentum space) operator V : 

A(V) = (x V£ + y V*) 5 LU , 
B(V) = (y-x) (V L A -V L B )5 LL ,, 
F(V) = y /xy(V*-V*) 5 LL , 

In case there is no off-diagonal disorder due to local lattice distortion because of size 
mismatch : 

S = ^^A(A- 1 )- 1/2 S RR/ A(A- 1 )- 1 / 2 r RR ,®X. 

R R! 
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The equation (8) is now exactly in the form in which recursion method may be applied. 
For ordinary recursion |1) is labeled by a particular L and has a representation as a 
1 x iV max column vector, with iV max being the size in augmented space. But in the case 
of block recursion, we deal with Lx V matrix and the length of |1) is L x V x iV max . We 
should note that at this point the above expression for the averaged <C Grl,ru{z) ^> 
or <C G LL i (k, z) 3> is exact. 

It is important to note that the operators A, B, F, J A , 3 B , J F , 6 A , 6 b and 6^ 
are all projection operators in real space (i.e unit operators in k- space) and acts on 
an augmented space basis only to change the configuration part (i.e. the cardinality 
sequence {C}). 

A||{C»=A 1 ||{C}> B\\{C}) = A 2 \\{C}) 8(R e {C}) F\\{C}} = A 3 \\{C ± R}) 
~3 A \\{C}) = J 1 \\{C}) J B \\{C}) = J 2 \\{C}) 5(R e {C}) ~3 f \\{C}) = J 3 \\{C ± R}) 
6 A \\{C}) = 0l \\{C}) 6 B \\{C}) = o 2 \\{C}) 5(R e {C}) 6 F \\{C}) = o 3 ||{C ± R}) 

The coefficients A\ — A%, J\ — J3 and o\ — 03 have been expressed in equation (9). The 
remaining operator S is off-diagonal in real space, but diagonal in k-space. 

In the real space representation: 

S \R,{C}) = J2 S( X ) \R + x,{C})- 

x 

Here x's are the nearest neighbour vectors. The operator S shifts the real-space site to its 
nearest neighbour position by the vector x an d configuration space remain unchanged. 
However, in the reciprocal space representation, it acts on an augmented space only to 
change the configuration part: 

s||{C}) = $>pH k -x) s( x ) ||{C- X }>. 

x 

Here the operator rigidly shifts the entire down-spin configuration by the vector \- The 
operation of the effective Hamiltonian is thus entirely in the configuration space. 

As long as we wish to obtain only the diagonal elements of the Green function, as is 
required for the local density of states or the spectral densities, the ordinary recursion 
as described by Haydock et al [1] suffices. However, for response functions we shall need 
the full Green matrix. This will be described in the following section. 

2.2. The configuration averaged current- current correlation function 

We shall quote here the result for the dominant term in the configuration averaged 
current- current correlation function reported by us earlier [5] . This correlation function 
is directly related to the optical conductivity. We note that the expression involves the 
full Green matrix in angular momentum space and not only its diagonal elements. 
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The expression for correlation function is 

d 3 k 



, z 2 ) > = / 

JBZ 



8tt 3 



Tr 



and the renormalized current term is given by : 



J eS (k, Zl ,z 2 ) = <j(k) >+2 



S(k,^ 2 )f(^)j( 1 )(k) + j«(k) f(^) E(k, 



+ S(k,z 2 ) f(^ 2 ) j (2) (k) f(*i) E(k,^ 



where 



f(z) = /^(z) = ( — 



C^__Cf _ /_L_ J_' 
At Af "Uf Af. 



J cff (k, zx, z 2 ) « G"(k, Zl ) > J cff (k, zx, z 2 ) f < G c (k, z 2 ) > 

(10) 



(11) 



and 



£ = 



G 



-i 



g is the virtual crystal Green function. The interested reader is referred to an earlier 
work [5] which derives these expressions in some detail. The main point in setting these 
equations out, is to note that in such calculations one needs the full Green matrix in 
angular momentum space. This is the main motivation for this work. 

2.3. Setting up the Block recursion 

The first step in setting up the Block recursion procedure is to systematically renumber 
the real-space basis with integers. An example on a square lattice is shown in figure 1. 

9 



10, 



11' 



12' 



13 

Figure 1. Systematic discrete numbering of nearest neighbour map on a square lattice. 



The nearest neighbour map is now generated by a systematic numbering of the states 
in augmented space as follows : 

(A) Real space formulation : We start with numbering \R, {0}) as 1, and then recursively 
generate the neighbours by acting on the states by S and F. Let us take an example of 
a square lattice : 
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(i) S acting on |1, {0}) = |1) gives four new neighbours |2, {0}) ... |5, {0}). The four 
real space neighbours of |1) are then |2), |3), |4) and |5). 

(ii) F acting on |1, {0}) = |1) is |1, {1}). This we number |6). 

(hi) S acting on |2, {0}) = |2) gives : |6,{0}), |7,{0}), |1,{0}) and |8,{0}). These we 

number |7), |8), |1) and |9). 
(iv) F acting on |2, {0}) = |2) is |2, {2}) = |10). 

We proceed exactly as above and finally obtain the nearest neighbour map matrix : the 
n-th column of whose m-th row gives the n-th neighbour of m. We show below the 
initial part of the nearest neighbour map for the above example : 

/ 2 3 4 5 6 \ 

7 8 1 9 10 

8 11 12 1 13 

v / 

The equivalences are : 
Augmented-space element => 1{0} 2{0} 3{0} 4{0} 5{0} 1{1} 6{0} 
Discrete numbering =^> 1 2 3 4 5 67 

Augmented-space element 7{0} 8{0} 2{2} 9{0} 1O{0} 3{3} 

Discrete numbering =^8 9 10 11 12 13 

(B) Reciprocal space formulation : 

In reciprocal space the procedure is even simpler, since the operators act only on the 
configuration part of the space. As before, we start with numbering \k, {0}) as 1 and 
then recursively generate the neighbours by acting on the states successively by S and 
F. Let us take the example of the square lattice : 

(i) S acting on |{0}) = |1) leaves it unchanged. The four neighbours of |1) are then 
|1), |1), |1) and |1). 

(ii) F acting on |{0}) = |1) is |{1}). This we number |2). 

(hi) S acting on |{1}) = |2) gives : |{2}), |{3}), |{4}) and |{5}). These we number |3), 

|4), |5) and |6). 
(iv) F acting on |{1}) = |2) is |{0}) = |1). 

We proceed as before and obtain the nearest neighbour map matrix. We again show 
below the initial part of the nearest neighbour map in reciprocal space : 

/ 1 1 1 1 2 \ 
3 4 5 6 1 
7 8 2 9 10 

v / 
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Cardinality sequence => {0} {1} {2} {3} {4} {5} {6} {7} {8} {1, 2} 
Discrete numbering =^123456789 10 

(C) The Block recursion : 

We now go over to a matrix basis of the form : {^jl^}, where J is the the discrete 
labeling of the augmented space states and L, L' labels the angular momenta. The inner 
product of such basis is defined by : 



J L" 

For a real space calculations on a lattice with Z nearest neighbours, we start the 
recursion with : 

*£L/ = Wg>, 5 JA + Wg, 5j, z+1 
While for a reciprocal space calculation we start with : 

where 

W « - 5 - 5 (12) 

™ll> - ~ TTi72 ° LL ' LL' - 7Tl/2 ° LL ' 



The remaining terms in the basis are recursively obtained from : 



L" J' L" L" 

®J,LL" B L"L> - 1^ 2^ H JL,J'L" ®J>,L"L> ~ 2^®J,LL" A L"L' ~ Z—i J,LL" B L"L' 
L" J' L" L" L" 

Orthogonalization of the basis gives : 

EEEE a©v *£W = E ^ 4t> (is) 

J L" J' L'" L" 

In matrix notation, where matrices are in angular momentum space : 

A (n) = ^nny 1 Q^HjJ' (14) 

J J' 

Next, we note that we had started with J max x L^ nax orthogonal basis set. The above 
procedure merely gives J max basis sets. We still have orthogonality conditions among 
the various columns of In order to impose these conditions consider : 

= EE^W^W - E^L-4t' - E*S5^< 

J' L" L" L" 



Block recursion and Green matrices 



9 



Construct L max column vectors (jy JL out of the L max columns of ^ j^ll 1 and set about to 
Gram- Schmidt orthonormalizing the set : 



LJ 

<p¥j = b 21 ^ j + b 22 ^j =► s 21 = e>SM 2 ]; ^ = EtfMi-^i 

LJ LJ 



n n—l 

= E *WS3 =► B nk = E^lVS (fc<n) ; ^ n = £>$MS " E S S* 

fc=l LJ LJ fc=l 

(15) 

We may now construct $j^ll> ou ^ °f V'jl an d note that -B n fe is indeed the matrix B^™ +1 ) 
we are looking for. 

The equations (14-15) show that we may calculate the matrices {A^,B^" +1 )} 
recursively, noting that = I and B^ ^ = 0. In this new basis, the Hamiltonian is 
block tri-diagonal and the Green matrix can be written as follows : 



G 



(n) 



E I - AW - B^+W G (n+ i) B(™ +1 ) 



i -i 



< G > = G(i) 



(16) 



3. Model Calculations 



(Aj Model on a square lattice : We shall first apply our methodology to a two-band 
model on a square lattice with the Hamiltonian : 



The concentration is taken to be x = 0.5. 

We have carried out block recursion for JV=10 levels. The termination was carried out 
as suggested earlier by Godin and Haydock [8] : we put {A< n \ B^+^j = {A^ N \ B^ N+1 ^} 
for n — N + 1, . . . , iV max and take G(jv max+ i) = (1/(E — iS)) I. In order to get a smooth 
density of states we had taken 5 = 0.01 and iV max = 10000. The elements of the 
Green matrix is shown in figure 2. The imaginary part of the diagonal elements give 
the projected density of states. Herglotz properties of the diagonal parts give rise to 
positive definite density of states. The off-diagonal part is relatively small and is not 
Herglotz. The projected density of states are symmetric, as is the imaginary part of the 
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Figure 2. The real and imaginary parts of the Green matrix for a 2 x 2 Hamiltonian 
model. The full lines refer to Gn, the dashed-dotted line to G22 and dashed lines to 

Gi2- 



off-diagonal element. The real parts of the matrix elements are also shown. These are 
related to the imaginary parts by the Kramers- Kronig relation. 

(B) The s-d model of a transition-noble metal alloy : Levin and Ehrenreich [20] and 
Gelatt and Ehrenreich [21] have introduced a simple two band model for transition-noble 
metal alloys. Physical effects like charge transfer between constituents will usually differ 
for the s-p conduction bands on one hand and the relatively narrow (i-bands on the other. 
Their model includes the conduction bands described together as a single band and the 
set of (i-bands also described by a single band and their hybridization. The model also 
takes into account the large widths of the conduction bands and the relatively narrower 
widths of the rf-bands. The following Hamiltonian has many (but not all) of the essential 
features : 

" = (18) 

The dominant disorder is taken to be in the terms e s and e^. The hybridization is taken 
between states at the same site only and the hopping terms are related by t s = atd- 
The sites R vary over the sites of a face-centered cubic lattice. As a model case we have 
taken the parameters shown in the following table : 



Constituent 


e s 


£d 


7 


a 


t d 


A 


1.5 


1.5 


0.2 


4.0 


-0.5 


B 


0.0 


0.0 


0.2 


4.0 


-0.5 



Table 1. Parameters for our calculation for A 50 B 5 o alloy on a fee lattice. 
(i) Calculations using real space block recursion : 
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Energy Energy 

Figure 3. The real and imaginary parts of the Green matrix for the s — d model on a 
fee lattice. The calculations were done by a real space block recursion. The full lines 
refer to Gu, the dashed-dotted line to G22 and dashed lines to G12. 

The figure 3 shows the real and imaginary parts of the Green matrix. These calculation 
are carried out through a real space block recursion technique. In general the qualitative 
features for the diagonal elements are similar to our square-lattice model. The main 
difference is that on a fee lattice the partial density of states, related to the imaginary 
part of the diagonal elements of the Green matrix, are no longer symmetric about the 
band centre. Consequently, the behaviour of the off-diagonal element is quite different. 
The total density of states is given by : 

n(E) = (1/tt) Im (G ss (E - i5) + 5G dd (E - i5)) . 

If, for example, the number of electrons in the constituents are 5 per atom per spin for 
A and 5.5 per atom per spin for B, the position of the Fermi energy is given by : 

/ dE n(E) = (n e ) = x n A + (1 - x) n B = 5.25 

J —00 

where n^ and n# are the number of valence electrons of A and B type of atom. 
(ii) Calculations using reciprocal space block recursion : 

We have carried out the block recursion in reciprocal space for the s — d model. The 
Green matrix in reciprocal space G(k, E) is the factor that arises in our earlier expression 
for the configuration averaged current-current correlation function. Its diagonal matrix 
element is related to the spectral functions for the different bands. In figure 4 we show 
the spectral functions along a given direction r — X in reciprocal space. The imaginary 
part of the off-diagonal matrix element is also shown in the figure. We note that the 
off-diagonal parts of the Green matrix has antisymmetric structure in its peaks, while 
the diagonal matrix elements are positive (as they represent spectral functions). 

Figure 5 shows the Green matrix elements calculated starting from the G(k, E) and 
carrying out a reciprocal space integration developed by us [23] as a generalization of 
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Figure 4. The matrix elements of G(k,E) on a fee lattice : the 11 clement (left), the 
22 element (middle), and the 12 element (right) along the r to X direction. 



the tetrahedron method proposed by Jepsen and Andersen [22] for crystalline systems. 
The close comparison between the figures 3 and 5 is a strong justification of the accuracy 
of the reciprocal space recursion. 



4. Optical conductivity of Ni 5 oPt 5 o alloy 

We shall now report a calculation of the optical conductivity of a 50-50 NiPt alloy. The 
initial electronic structure calculations require input of the potential parameters from 
the pure Ni and Pt. This we have taken from a TB-LMTO calculation on the pure 
materials. We use these as an input for a self-consistent calculation of the electronic 
structure using our LDA-self consistent Augmented space recursion technique (ASR). 
We first calculate the Green matrix in reciprocal space. The results are shown in figure 6. 
The imaginary part of the diagonal elements yield the spectral functions. The figure 6 




-10 o 
Energy 



10 
Energy 



Figure 5. The real and imaginary parts of the Green matrix for the s — d model 
on a fee lattice. The calculations are done by a k-space block recursion followed by a 
reciprocal space integration [23]. The full lines refer to Gn, the dashed-dotted line to 
G22 and dashed lines to G12. 
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Figure 6. The imaginary and real parts of (l/9)Tr Gll{^, E) (left column) and the 
imaginary and real parts of Gll'(E) (right column) for Ni 5 oPt 5 o alloy. 



also shows the off-diagonal matrix elements. Most of these elements are very small as 
compared with the spectral functions and are not positive definite as expected. The 
density of states for the alloy is shown in figure 7. We obtain the density of states by 
integrating the spectral function over the Brillouin zone, using the tetrahedron method 
generalized for disordered alloys by us [23]. The density of states for NiPt is particularly 
difficult to reproduce accurately by real space recursion, because of the sharp structure 
straddling the Fermi energy. However, the calculation in reciprocal space followed by 
Brilluoin zone integration reproduces the density of states in very good agreement with 
earlier work using KKR-CPA [24]. 

The calculation of the optical conductivity involves the calculation of the four current 
terms } AA , j AB , j BA and j BB . Ideally one should calculate these out of two atoms of 
the type AA, AB, BA or BB embedded in the disordered alloy. However, as a first 
approximation we have obtained these from the pure metals and the ordered alloy. 
Having obtained the current terms we use the scattering methodology described in [5] 
to obtain the correlation function. 



In figure 8 we show the result for the correlation function (full lines) for a 50-50 NiPt 
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Energy (eV) 



Figure 7. The density of states for a 50-50 NiPt alloy. The dashed line shows the 
Fermi energy. 

alloy : 

S(u>) = 1 J dE, [S(Ef, Ef) + S(E^, £ 2 -) - S(E+, E%) - S(E£, E+) _ . 

where E 2 = E 1 + w, E ± = E ± i0 + . For comparison we have also shown (as dashed 
lines) the joint density of states for the same alloy 

J(u) = J dE n v {E) n c {E + uj). 

The former calculation requires the full Green matrices as well as the full matrix self- 
energies (matrices in the angular momentum space). We have carried out the calcula- 
tions of these matrices via our block recursion technique. The joint density of states is 
modulated by the current terms. The main structures of the joint density of states are 
retained in the correlation function, but the relative heights of the principal peaks are 
modified. 



5. Conclusion 

In this communication we have described a block recursion in augmented space suitable 
for calculations of the Green matrices. The recursion is set up both with real space 
augmented by the configuration space of the alloy or the reciprocal space augmented 
with the configuration space. For the latter case we have coupled it with a Brillouin 
zone integration scheme which is a generalization of the tetrahedron method developed 
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Figure 8. The joint density of states J(w) (dashed line) and the correlation function 
S(u) (solid line) for 50-50 NiPt alloys. 

earlier for crystalline systems. The Green matrices are essential for the calculation of 
the response functions and effective current terms which are related to the self-energy 
matrices. We propose to use these techniques in our future applications. 

References 

[1] R. Haydock, V. Heine and M. J. Kelly J. Phys. C: Solid State Phys. 5 2845 (1972). 

[2] T. Saha, I. Dasgupta and A. Mookerjee, J. Phys.: Condens. Matter 6 L245 (1994). 

[3] T. Saha, I. Dasgupta and A. Mookerjee, J. Phys.: Condens. Matter 8 1979 (1996). 

[4] A function of a complex variable f(z) is called Herglotz if : all singularities of f(z) lie on the real z 
axis and Im f(z) = — sgn(Im z) for all z off the real axis. 

[5] K. K. Saha and A. Mookerjee, e-print cond-mat/0403456 and Phys. Rev. B 70 (2004) (in press). 

[6] A. Alam and A. Mookerjee, e-print cond-mat/0407512. 

[7] T. J. Godin and R. Haydock, Phys. Rev. B 38 5237 (1988). 

[8] T. J. Godin and R. Haydock, Phys. Rev. B 46 1528 (1992). 

[9] A. Mookerjee, J. Phys. C: Solid State Phys. 6 1340 (1973). 
[10] A. Mookerjee and R. Prasad, Phys. Rev. B 48 17724 (1993). 
[11] I. Dasgupta, T. Saha and A. Mookerjee, Phys. Rev. B 51 3413 (1995). 
[12] T. Saha and A. Mookerjee, J. Phys.: Condens. Matter 8 2915 (1996). 
[13] A. Chakrabarti and A. Mookerjee, J. Phys.: Condens. Matter 14 3211 (2002). 
[14] K. K. Saha and A. Mookerjee, J. Phys.: Condens. Matter 16 1409 (2004). 
[15] S. Ghosh, N. Das and A. Mookerjee, Int. J. Mod. Phys. B 21 723 (1999). 

[16] P. Biswas, B. Sanyal, M. Fakhruddin, A. Haider, A. Mookerjee and M. Ahmed, J. Phys.: Condens. 
Matter 7 8569 (1995). 

[17] P. Biswas, B. Sanyal, A. Mookerjee, A. Huda, N. Chowdhury, M. Ahmed and A. Haider, Int. J. 

Mod. Phys. B 11 3703 (1997). 
[18] O. K. Andersen, Phys. Rev. B 12 3060 (1975). 
[19] O. K. Andersen and O. Jepsen, Phys. Rev. Lett. 53 2571 (1984). 



Block recursion and Green matrices 



16 



[20] K. Levin and H. Ehrenreich, Phys. Rev. B 3 4172 (1971). 
[21] C. D. Gelatt, Jr. and H. Ehrenreich, Phys. Rev. B 10 398 (1974). 
[22] O. Jepsen and O. K. Andersen, Solid State Commun. 9 1763 (1971). 
[23] K. K. Saha, A. Mookerjee and O. Jepsen, c-print cond-mat/0405175. 

[24] F. J. Pinski, Ginatempo, D. D. Johnson, J. B. Staunton, G. M. Stocks and B. L. Gyorffy, Phys. 
Rev. Lett. 66, 766 (1991). 



